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Objective: Extend a 3-D combustor program: 

• To predict pollutant emissions of soot and N0 X 

• To include the influence of soot, CC^, and on radiation heat 
transfer 

• To extend the two-step hydrocarbon oxidation mechanism to a more 
detailed four-step scheme 

Background : EPA regulations for aircraft gas turbines make it imperative that 

soot and NO emissions be controlled. The use of alternative synthetic 

A 

fuels, recently being considered, results in significant increases in 
soot production. Thus, an improved understanding of the physical and 
chemical processes governing soot production is necessary. Since these 
processes are not well understood, only an approximate modeling of soot 
production is currently possible. Lack of general 3-D computer codes has 
also hindered research in this area. In the present program, soot and 
NO emissions were computed by using an existing general 3-D computer 
program. This will significantly aid in the development of soot and 
NO production models in practical systems. 

Approach: A general 3-D combustor performance program developed by Garrett 

was extended to predict soot and N0 X emissions. The soot formation and 
oxidation rates were computed by quasi-global models, taking into account 
the influence of turbulence. Radiation heat transfer was computed by the 
six-flux radiation mode. The radiation properties include the influence 
of C0 2 and H^O in addition to soot. N0 X emissions were computed 
from- a global four-step hydrocarbon oxidation scheme and a set of 
rate-controlled reactions involving radicals and nitrogen oxides. 
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Results: Computations performed for a plug flow reactor show the four-step 

scheme to be far superior to the two-step scheme in predicting 
temperature and species concentrations. 

• Computations were performed for idle, cruise, and takeoff conditions of 
a JT8D combustor. These showed that the present model is capable of 
producing reasonable predictions of smoke and NO emissions and of 

A 

the wall radiation flux. 
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COMPUTATIONS OF EMISSIONS USING 
A 3-D COMBUSTOR PROGRAM 


OBJECTIVE: 

EXTEND AN EXISTING 3-D COMBUSTOR PROGRAM: 

• TO PREDICT POLLUTANT EMISSIONS OF SMOKE AND NO* 

• TO INCLUDE THE INFLUENCE OF SOOT. C0 2 . AND H 2 0 ON 
RADIATION HEAT TRANSFER; AND 

• TO EXTEND THE TWO-STEP HYOROCARBON OXIDATION 
MECHANISM TO A MORE DETAILED FOUR-STEP SCHEME 


APPROACH 

PROGRAM INVOLVED FOUR TASKS 

• TASK I - FORMULATION OF THE METHOD 

• TASK II - COMPUTER CODING 

• TASK III - COMPUTATION OF TEST CASES 

■ IDLE. CRUISE. AND TAKE-OFF CONDITIONS FOR A 
JTBO COMBUSTOR 

• TASK IV - REPORTING AND DOCUMENTATION 
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EXISTING 3-D COMBUSTOR PROGRAM 


• GENERAL PROGRAM FOR 

■ RECIRCULATING. SWIRLING. TURBULENT. REACTING FLOW 

■ GASEOUS AND/OR LIQUID FUEL COMBUSTION 
> RADIATION NEAT TRANSFER 

• GENERAL TRANSPORT EQUATION 

div fc>u£ - GRAD <t > ) * S <t> 

• VARIABLES: UT>. k. r. I. M| U . Mco (Mo,. Mflb. Mr? 01- 
h |T. p\. R x . R y . R z SPRAY DYNAMICS/COMBUSTION 

• PHYSICAL MODELS: 

■ TURBULENCE: k-t MODEL 

■ CHEMISTRY: 2-STEP REACTION SCHEME 

• CHEMICAL REACTION RATE: MODIFIED EBU MODEL 

■ RADIATION: SIX-FLUX MODEL 

• NUMERICAL 

• FINITE-DIFFERENCE ITERATIVE METHOD 

■ SUITABLE FOR COMPLEX GEOMETRIES 

■ SUITABLE FOR NON-UNIFORM GRID SPACING 


FOUR-STEP HYDROCARBON OXIDATION MECHANISM 


• TWO-STEP SCHEME IN ORIGINAL 30 PROGRAM: 

CxHy + O2 - CO + H2O 

CO + O2 — CO2 

• FOUR-STEP SCHEME: 

C X Hy - C X Hy _2 + H2 
CvHw? + 0 2 - CO + H2 
CO + 0 2 - CO2 

H2 + O2 — H2O 

• FOUR-STEP SCHEME DESCRIBES THE FOLLOWING BASIC STEPS 
OF HYOROCARBON OXIDATION 

■ TRANSFORMATION OF HYOROCARBON FUEL INTO INTERMEDIATE 
HYDROCARBONS AND HYDROGEN WITH LITTLE RELEASE OF ENERGY 

■ OXIDATION OF INTERMEDIATES TO CO AND H2 

■ OXIDATION OF CO TO CO2 

■ OXIDATION OF H2 TO H2O 

• TWO ADDITIONAL EQUATIONS FOR: C x H y .2 AND H* 
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N0 X EMISSIONS 

• SPECIES CONSIDERED: C x H y . C x H y . 2 . CO. N. NO. NO 2 . H. 

H 2 . OH. 0. C0 2 . H 2 O. O 2 . N 2 

• TRANSPORT EQUATIONS FOR: C x H y . CO. C x H y . 2 , N. NO. NO 2 . 
H. H 2 . OH, 0 

• C0 2 . H 2 0. O 2 FROM C. H. 0 ELEMENT CONSERVATION 

• Mn2 = 1-IM 

• SOURCES OF SPECIES COMPUTED FROM MODIFIED EBU MOOEL 

• C x H y . C x H y . 2 . H 2 . AND CO SOURCES FROM 4-STEP 
REACTION SCHEME 


N0 X EMISSIONS 


• SOURCES OF N. NO. NO 2 . H. H 2 . OH. 0 FROM REACTION MECHANISM: 


0 + N 2 = N + NO 
NO + 0 = O 2 + N 
N + OH = NO + H 
N 2 + 0 2 = N + NO 2 
N 2 + 0 2 = NO + NO 
NO + NO = N + NO 2 
M + N0 = M + N + 0 
M + N0 2 = M + 0 + NO 
M + N0 2 = 0 2 + N + M 
NO + 0 2 = NO 2 + 0 


H + N0 2 = NO + OH 
H 2 + 0 2 = OH + OH 
OH + H 2 = H 2 0 + H 
OH + OH = 0 + H 2 0 
0 + H 2 = H + OH 
H + 0 2 = 0 + OH 
M + 0 + H = 0H + M 

M + 0 + 0 = 0 2 + M 

M + H + H = H 2 + M 

M + H + OH = H 2 0 + M 


STRONGLY COUPLED NONLINEAR STIFF EQUATIONS 

PRATT'S CREK PROGRAM FOR STIFF EQUATIONS AOAPTEO TO 
3-0 COMBUSTOR PROGRAM 
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SOOT EMISSIONS 


• VARIOUS STEPS INVOLVED IN SOOT FORMATION/OXIDATION NOT QUANTIFIED 

• HENCE QUASI-GLOBAL MODELS CONSISTING OF THREE STEPS 

■ SOOT NUCLEI FORMATION 

■ SOOT PARTICLE FORMATION 

■ SOOT PARTICLE OXIDATION 

• TRANSPORT EQUATIONS FOR NUCLEI ANO SOOT CONCENTRATIONS WITH 
SOURCES CONTAINING REACTION RATES 

• REACTION RATES CONTAIN FORMATION AND OXIDATION RATES DEPENDENT ON: 

■ UNBURNT FUEL AND 0 2 CONCENTRATION 

■ TEMPERATURE 

■ TURBULENCE PARAMETERS d k (SIMILAR TO EDDY-BREAK-UP MODEL) 

• TWO SIZES OF SOOT PARTICLES CONSIDERED: ONE EQUATION FOR EACH SIZE 


RADIATION HEAT TRANSFER 

• SIX-FLUX RADIATION MOOEL - SCHUSTER ANO HAMAKER 

• MAJOR CONTRIBUTORS TO RAOIATION = SOOT. C0 2 . H 2 0 

• RAOIATION PROPERTIES DEPENDENT ON SOOT. C0 2 . H 2 0 CONCENTRATIONS 
? TOTAL PROPERTIES (INTEGRATED OVER WAVELENGTH| 

• C0 2 - H 2 0 EMISSIVITY CORRELATIONS INCLUDING OVERLAP ANO 
PRESSURE CORRECTION FACTORS 

*C+W = *C + *W C W — C W Af 
»C+W = Ln M- r C+Wl 1 s 

• SOOT RADIATION PROPERTIES 

■ ABSORPTION COEFFICIENT COMPUTEO FOR BOTH PARTICLE SIZES 
AND ADDED TOGETHER 

• GAS-SOOT MIXTURE, 

> ABSORPTION COEFFICIENTS OF GAS AND SOOT ADDED 
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RESULTS 


• FOUR-STEP SCHEME 

■ PLUG FLOW REACTOR FOR LEAN. STOICHIOMETRIC. AND 
RICH PROPANE FLAMES 

• CONSIDERABLE DEVIATIONS OF TWO-STEP PREDICTIONS 
FROM MEASUREMENTS 

• FOUR-STEP PREDICTIONS AGREE CLOSELY WITH 
MEASUREMENTS 

■ SOME DISCREPANCY IN FOUR-STEP H Z PREDICTIONS AT 
HIGHER EQUIVALENCE RATIOS 


LEAN C 3 H 3 FLAME (0 = 0.12) 
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TEMPERATURE. K 


LEAN C3H8 FLAME (0 




AXIAL DISTANCE. CM 


AXIAL DISTANCE. CM 


LEAN C 3 H 8 FLAME (0 = 0 . 12 ) 


V 

\ CH INTERMEDIATES VN 


AXIAL DISTANCE. CM 



AXIAL DISTANCE. CM 







JT8D-17 COMBUSTOR 












PREDICTED EMISSIONS INDEX WITH FOUR-STEP 
HYDROCARBON OXIDATION SCHEME 



EMISSION INDEX 
Gm OF EMISSIONS/Kg OF FUEL 


CONDITION 

SMOKE 

N0 X 

IDLE 

0.056 |0.6) 

M3 

CRUISE 

1.3 

13 

TAKEOFF 

1 2 (2.8) 

27 (24.4) 


NOTE: VALUES IN PARENTHESES ARE EXPERIMENTAL MEASUREMENTS 

• PREDICTED RADIATION FLUXES IN THE SAME RANGE AS 

MEASUREMENTS FOR SLIGHTLY DIFFERENT OPERATING CONDITIONS 


CONCLUSIONS 


• MODEL PRODUCES REASONABLE RESULTS FOR EMISSIONS 
AND RADIATION FLUX 


• LACK OF EXPERIMENTAL DATA PRECLUDES MORE DETAILEO 
MODEL VALIDATION ' 
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